Luttinger liquid superlattices 
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We calculate the correlation functions and the DC conductivity of Luttinger liquid superlattices, 
modeled by a repeated pattern of interacting and free Luttinger liquids. In a specific realization, 
where the interacting subsystem is a Hubbard chain, the system exhibits a rich phase diagram with 
four different phases: two metals and two compressible insulators. In general, we find that the 
effective low energy description amalgamates features of both types of liquids in proportion to their 
spatial extent, suggesting the interesting possibility of 'engineered' Luttinger liquids. 
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In recent years, new experimental techniques have 
made it possible to grow nanostructures which are topo- 
logically one-dimensional, such as quantum wires and 
Carbon nanotubes. However, care must be taken when 
invoking existing models to discuss their electronic prop- 
erties, since inhomogeneities must be taken into account 
in a fundamental way. Consider, for instance, the Lut- 
tinger liquid (LL), which is the standard model for low- 
energy phenomena involving interacting electrons in one 
dimension [pj . The absence of conductance renormaliza- 
tion in long high-mobility GaAs wires g has been ex- 
plained in terms of a usual LL (representing the wire) 
in contact with a non-interacting LL at each of its ends 
(representing the Fermi liquid leads); that is, overall, the 
system can be thought of as an inhomogeneous Luttinger 
Liquid (ILL). LL's with different inhomogeneity profiles 
have also been used in the context of the fractional quan- 
tum Hall effect (FQHE), to describe transitions between 
edge states at different fillings (|fl, or between an edge 
state connected to a Fermi liquid B . A different class of 
inhomogeneous systems is represented by superlattices 
(SL's) and multilayers. By varying the relative thick- 
nesses of the repeating unit in magnetic metallic multi- 
layers, fascinating properties such as exchange oscillation 
and giant magneto-resistance (GMR) have been found 
pL The interplay between electron correlations and a 
SL structure (or layering) can therefore lead to collective 
properties quite distinct from those of each of its con- 
stituents. Furthermore, the ability to manipulate physi- 
cal properties by choosing an appropriate spatial modu- 
lation opens the way for a whole new set of 'engineered' 
materials. 

With this in mind, our purpose here is to discuss the 
properties of a one-dimensional SL made up of a peri- 
odic arrangement of two long and perfectly connected 
LL's, one interacting and the other free. Accordingly, 
the low-encrgy properties of this Luttinger liquid super- 
lattice (LLSL) are described by generalizing the usual 
bosonized Hamiltonian fplnMI as follows 
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where the sum extends over separated charge- (v = p) 
and spin- (y = a) degrees of freedom, each of which with 
layer-dependent parameters u u {x) and K v {x); these de- 
termine, respectively, the velocity of elementary excita- 
tions and the algebraic decay of correlations functions for 
each degree of freedom. For x on the free layer one has 
K v (x) = 1 and u u (x) = vf, the Fermi velocity, whereas 
for x on the repulsive layer K v {x) and u v {x) become the 
usual uniform LL parameters. For definiteness, we will 
often speak of a 'Hubbard superlattice' (HSL), where the 
interacting layer is taken as a Hubbard model with hop- 
ping t and on-site repulsion U [PlJlOfl; a weak coupling 
perturbation theory similar to that of the homogeneous 
model can be used to show that Eq. (Q) indeed describes 
the low energy and small momentum sector of the dis- 
crete model with long layers 1 1 1 1 . In the homogeneous 
case, the dependence of the LL parameters on both the 
density and U has been determined by recourse to the 
exact solution ]l^-fiq]. With respect to magnetic prop- 
erties, the SL structure does not break spin SU(2) sym- 
metry, so that the inhomogeneous K a is still expected 
to renormalizc to K* = 1. The spin sector stiffness is 
therefore unrenormalized as in the homogeneous system 

©■ 

The boson phase fieids <& v are related to the charge and 
spin densities, p and er, through \/2d x ^ > ^(x)/TT = v, while 
0;, is such that <9 æ 0„ is the momentum field conjugate 
to $,,: [<& v (x), dy@ u '(y)} = i5yy5(x-y). § u and <d u are 
dual fieids, since they satisfy both 



dt$ v = u v (x)K v (x)d x Q v 



(2) 



and the equation obtaincd through the replacements 
®v — * ©i/j 9^ — ► <ÉV, and K v — » 1/K V . These equa- 
tions can be uncoupled to yield 



du*v - u v K„d x ( j^d x $„ ) = o, 



(3) 
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and a dual equation for &„. 



The equations of motion are subject to the continuity 
of $i, and 0„ [jl^,[l^] (which ensures the continuity of 
the fermionic field). Since their time derivatives are also 
continuous, Eq. (0) and its dual yield, as additional con- 
ditions, the continuity of (u v j K v ) d x <& v and u v K v d x Q v 
at the contacts. Thus, both charge and spin currents 
j u = V2d t & u /n are conserved, since we neglect Umk- 
lapp processes and spin backscattering. 

The Hamiltonian p) is straightforwardly diagonalized 
by a normal mode expansion 



Q v {x,t) = i 22 
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where b' pv are boson creation operators (p > 0). (f>o tU (x) 
and 6q )V (x) are the zero mode fimetions which, in the 
homogeneous case, are given by (f>o tV (x) = N v ^- , 
6*0,1/ (x) — Jv^-t where N v and J v are the (charge and 
spin) number and current operators. Besides, in this 
case 7„ = ttu v K v J,j/ L and t v — n(u u /K v )N u /L. How- 
ever, in a LLSL the inhomogeneity will induce a modu- 
lation of the charge (but not of the spin) density of the 
system. Thus, one needs to introduce in general layer- 
specific number and current operators. Since each layer 
is a LL, the variations across it are A</>o.„ = ttN\ v and 
&-6q,v = kJ\v, where Å = or U, depending on whether 
it is a free or interacting layer, respectively. <j>o,v{x) and 
6*0,1/ (x) will then be linear continuous functions of x, with 
slopes given by the layer number and current operators 
(we omit the expressions for brevity). Analogously, from 
the equations of motion (g), we obtain "f\ y and t\ v . 

In order to find the equilibrium value of the density in 
each layer, one needs to equate their chemical potentials 



Mofao) = Mo (n + £(n - nu)) = nu(nu), 



(G) 



where n = N/L is the total electron density, £ — Ljj/Lq 
and u\ and n\ are the chemical potential and density of 
each layer, respectively. For definiteness, we have de- 
termined the charge profile in a HSL using the exact 
expression for (J,u{nu) [Q- We found that the charge 
tends to accumulate in the free layer. This is rather in- 
tuitive, since electrons decrease their mutual repulsion 
energy by flowing into the free layer. This was observed 
in numerical studies of the HSL JTn] . Of course, such a 
charge inhomogeneity will be strongly suppressed with 
the inclusion of long-range Coulomb interactions, which 
are absent in a Hubbard model description. 

The HSL has a very rich phase diagram. For n < 1 
the system is always metallic. For n > 1, however, we 
observe four different phases, two metallic and two insu- 
lating, each characterized by its charge profile, as shown 
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FIG. 1. Phase diagram of a Hubbard superlattice show- 
ing two metallic (M) and two insulating (I) phases [i = 1) 
(U c /t « 3.2309 and n c = (2 + £)/{l + tj). 

in Fig. P. The two insulating phases correspond to ei- 
ther nu = 1 (Mott insulator) or n$ = 2 (band insulator). 
Note that, in each case, one type of layer is by itself in- 
sulating while the other is metallic, the overall insulating 
character being a consequence of the ID structure ('re- 
sistors in series'). Therefore, both insulating phases are 
gapless, except at the phase boundary indicated by the 
dashed line in Fig. ni where both nu — 1 and n$ = 2 and 
the system exhibits a Mott-Hubbard gap. The density at 
this line is thus n c = (2 + £)/(l +£) @. The two metal- 
lic phases differ in the density of the interacting layer: 
the chemical potential can fall in either its upper or its 
lower Hubbard band. For values of U > U c æ 3.2309£, 
the system cannot sustain the upper metallic phase. This 
value is independent of the 'aspect ratio' £. Note that, 
for U < U c , the HSL is always gapless. 

While the topology of the phase diagram of Fig. [l] is 
specific to a HSL, we expect that several of its features 
should be generic to other LLSL's. In particular, the 
'division of labor' between the two subsystems, where 
one is responsible for the insulating behavior whereas the 
other renders the system compressible, is reflected in the 
weighted form of the SL compressibility 
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where nu = ^Kp/nup and Ko = 2/ttvf are the compress- 
ibilitics of the interacting and free layers respectively. 

The SL structure also affects the velocity of excitations. 
Forp <§C tt / (Lu + Lq) , the dispersion relation of the LLSL 
is linear, with effective velocities 
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where A^ = K v + K^ 1 . Clear ly, c v ^u„ as £ — > oo, and 
Cv —* Vf as £ — > 0. Furthermore, as one approaches the 
insulating phase from the low-density region (see Fig. [l]), 
c p — > as a result of u p — > in the interacting layer. As 
in the homogeneous system, the LL description breaks 
down whenever a gap opcns in the charge or spin sector 
of either layer. In the HSL case, this happens in both 
insulating regions of Fig. pi. Note, however, that the de- 
termination of the phases through Eq. (ph does not rely 
on the LL description. 

We now focus on the correlations. The T = asymp- 
totic behavior (i.e., for well separated x and y) of charge 
and spin correlations is given by 
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where (j)^{x) ~ 


k F x — 4>o^p(x) and the LLSL effective ex- 


ponent is 
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here a v is a function of system parameters and the layer. 
Similarly, correlation functions for singlet and triplet su- 
perconducting pairing are 
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where K p — f(l/K p ). One should note that the corre- 
lation functions depend not only on the difference x — y, 
but also on the actual positions x and y, through the 
zero mode functions. Their effect will be to generate the 
usual spatial oscillations present in homogeneous LL's. 
However, due to the inhomogeneous density profile, their 
period will vary from layer to layer, reflecting the layer- 
dependent Fermi wave-vectors; this is akin to the os- 
cillatory behavior of the exchange coupling in magnetic 
metallic multilayers p|L0||. In spite of the presence of 
effective exponents K* and K p , the condition for domi- 
nant superconducting correlations reduces to the one for 
homogeneous systems, namely K p > 1. The dominant 
term in the charge and spin correlation functions is the 
second one, which in the homogeneous case corresponds 
to the 2k F contribution. This predominance, however, 
may be superseded by the behavior of the amplitude A\ , 
as discussed in Ref . Q . 

In Fig. H, the correlation exponent K* of a HSL is 
shown as a function of Alling, for £ — 1. Both metallic 
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FIG. 2. The correlation exponent K* of a Hubbard super- 
lattice as a function of the total electron density n for £ — 1 
and different values of U. 

phases are characterized by 1/2 < K* < 1. On the low 
density side, K* approaches a value larger than 1/2 as 
n — ► 0, which depends on £ but not on U . This is a 
feature unique to the LLSL. Note from Eq. ( pd| ) that K* 
interpolates monotonically between 1 (the free value) and 
K p (the interacting layer exponent) as £ is varied from 
to co. This illustrates a general feature of the LLSL, 
namely: by varying the 'aspect ratio' £, one can fine-tune 
a physical property to a specified value. 

Finally, we discuss the transport properties of a LLSL. 
In the presence of an applicd electric field the equation 
of motion for $ p becomes |l7| , [18| 
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$ p (x,t) =-eE(x,t). (13) 



The nonlocal conductivity is given by 

2g a 



cr(x,y,t) 



d t G(x,y,t) 



(14) 



where g — e 2 /h is the conductance quantum and 
G(x,y,t) = -i 0(t) ([$p(x,t),$ p (y, 0)]) is the bosonic 
Green's function. 

We first calculate the Drude conductivity, which 
gives the current response to a uniform electric field: 
lim^^o a(q — 0, u>) pl| . A straightforward calculation 
yields 



a(q = 0,w) = 2g c p K* 6(u), 



(15) 



The delta function coefficient is the Drude weight. It has 
the same form as for the homogeneous case |n2| , but with 
the effective velocity and effective exponent replacing the 
corresponding uniform quantities u p and K p . By plug- 
ging in the results from Eqs. (0) and (|il[), one recognizes 
the conductivity of resistors connected in series. 
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FIG. 3. The Drude weight of a Hubbard superlattice as a 
function of n for £ = 1 and different values of U. 



engineered to suit a particular purpose, in a way remi- 
niscent of modulation-doped semiconductor heterostruc- 
tures and magnetic multilayers. Whether this will prove 
feasible, however, remains to be seen. 

In summary, we have considered Luttinger liquid su- 
perlattices made up of a periodic arrangement of frcc and 
rcpulsivc Luttinger liquids. Due to the space-dependent 
properties of the system, a non-homogeneous charge pro- 
file ensues. A specific realization of such a system, a 
Hubbard superlattice, was investigated in detail and its 
phase diagram was shown to exhibit two metallic phases 
and two peculiar compressible insulating ones. 
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We have plotted the Drude weight of a HSL as a func- 
tion of n for l = 1 and several values of U in Fig. f| The 
plot shows the re-entrant behavior as a function of n for 
U < U c . Furthermore, the Drude weight dips to zero 
upon approaching the insulating regions as a result of 
the vanishing charge velocities u p — » (Mott insulator) 
and vf — ► (band insulator). 

A more common experimental situation occurs when a 
field is applied to a finite region of the sample. In this 
case, the inverse order of limits applies lim g ^o 0"(<Z, w — 
0) , and we obtain the Landauer conductance J2M . In the 
LLSL 



a(q,0J = 0) = 2g o K;S(q). 



(16) 



Once again, the result is in close analogy with the ho- 
mogeneous system [|22|, where the SL interaction expo- 
nent has replaced the homogeneous one. Both Drude 
and Landauer responses, therefore, can be modulated 
by changing the SL spacing. However, the interaction 
renormalization of Eq. (Uq) is not revealed in usual DC 
conductance measurements, where it is masked by the 
presence of the Fermi liquid leads JL7J. In a LLSL of 
length L, only by going to frequencies of the order of the 
inverse traversal time u> > u p /L can the influence of the 
K* exponent be felt g. 

In closing, we would like to highlight the faet that 
the low energy, long wavelength properties of a LLSL 
can be, in effect, subsumed into a few effective parame- 
ters, in close analogy with the usual homogeneous LL dc- 
scription. These effective parameters, on the other hånd, 
turn out to be weighted averages of the underlying sub- 
system properties, in rough proportion to their spatial 
extent [see, e.g., Eqs. (^), (||) and Æj])]. Such a 'tem- 
pered Luttinger liquid' description suggests the interest- 
ing possibility of creating SL structures with properties 
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